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Possible Odd-Frequency Superconductivity in Strong- Coupling Electron-Phonon Systems 
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A possibility of the odd-frequency pairing in the strong-coupling electron-phonon systems is discussed. 
Using the Holstein-Hubbard model, we demonstrate that the anomalously soft Einstein mode with the 
frequency a>E ci>c {f^c is the order of the renormalized bandwidth) mediates the i-wave odd-frequency 
triplet pairing against the ordinary even-frequency singlet pairing. It is necessary for the emergence of the 
odd-frequency pairing that the pairing interaction is strongly retarded as well as the strong coupling, since 
the pairing interaction for the odd-frequency pairing is effective only in the diagonal scattering channel, 
(w„, —cj„) — » {(jJ,,', with ai,,/ = <jj„ > aj^- Namely, the odd-frequency superconductivity is realized in 
the opposite limit of the original BCS theory. The Ginzburg-Landau analysis in the strong-coupling region 
shows that the specific-heat discontinuity and the slope of the temperature dependence of the superfluid 
density can be quite small as compared with the BCS values, depending on the ratio of the transition 
temperature and ojc- 
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1. Introduction 

A novel symmetry of superconducting pairing, which is 
odd in imaginary time of the relative coordinate of the Cooper 
pair, has been investigated for several decades. The so-called 
odd-frequency (OF) superconductivity was first proposed by 
Berezinskii as a possible triplet pairing in the A-phase of su- 
perfluid ^He.'' Two decades later, Balatsky and Abrahams re- 
visited the singlet type of the OF pairing in the context of 
cuprates superconductors.^^* Since then, the OF pairing has 
been discussed in a wide variety of theoretical models, e.g., 
the Kondo lattice,^"^' the square-lattice'"' and the triangular- 
lattice""'-'' Hubbard, and r-7'"'' models. A universal feature 
of superconductivity so far stimulates to discuss the OF pair- 
ing in connection with experimental reality in the heavy- 



fermion systems 



15, 16) 



superconducting junctions 



17-30) 



vortex 



core,^''^^' proximity effect of superfluid ^He^^' and the cold 
atoms.-'"'' 

Aside from the intensive studies, the homogenous state of 
the OF pairing has struggled against the fundamental difficul- 
ties such as the thermodynamic instability^' ^^•■'*' and the un- 
physical Meissner effect.''' These difficulties have hampered 
further investigation on the property of the superconducting 
state except a determination of the transition temperature Tc- 
Recently, it is pointed out that these difficulties can be settled 
by the appropriate treatment of the gap function with retar- 
dation in the path-integral framework. ^^••'^' As a result, the 
property of the OF state is ready to be examined without any 
theoretical ambiguity. 

In the OF pairing, a strong on-site Coulomb repulsion is 
avoided by the off-time pairing with a vanishing equal-time 
pair amplitude. In order to gain an attraction efficiently and 
overcome an ordinary even-frequency (EF) pairing, a strongly 
retarded interaction, i.e., with strong frequency dependence, 
as well as a repulsive part of the interaction varying slowly 
in frequency may be preferable. Following this scenario, the 



realization of the OF pairing was examined by the mecha- 
nism of exchanging the critical antiferromagnetic fluctuation 
near the magnetic quantum critical point and the massless spin 
wave in the antiferromagnetic phase. '^* 

In this paper, we focus on peculiar electron-phonon sys- 
tems such as Skutterudite''^' and yS-pyrochlore"*"' families as 
another possible candidate for the OF pairing. These sys- 
tems are characterized by a cage-like structure with a large 
anharmonicity of an ionic potential where relatively isolated 
phonons constitute well-defined Einstein modes.""' They are 
expected to be particularly robust against strong electron- 
phonon coupling inherent from the cage-like structure. The 
investigations based on the Holstein-Hubbard model for such 
systems concluded that dynamical effect leads to the further 
softening of the Einstein mode and the conduction band is 
heavily renormalized near the half filling. ''^^'*' From these 
conclusions, the effective attractive interaction works in a 
very narrow range e < and the strength of the pair- 
ing attraction falls into the strong-coupling region. Namely, 
/lep = pfVo ~ ,?^/weWc > 1, where g, we and coc are 
the renormalized electron-phonon coupling, the Einstein fre- 
quency and the bandwidth, respectively, with the condition 
We ^ <^c- It is the situation that is quite suitable for the OF 
pairing, although only the possibility of the EF pairing has 
been so far investigated extensively.''-'"^^' 

The purpose of this paper is two folds: (i) to investigate the 
favorable conditions for the OF pairing in the strong-coupling 
electron-phonon systems, (ii) to elucidate the properties of the 
s-wave triplet as the simplest case of the OF pairing. In §2, we 
give the Holstein-Hubbard model and the formulation for the 
OF pairing state. The discussion of the structure of the pairing 
interaction and the properties of the OF phase are given in §3. 
The last section is devoted to the Ginzburg-Landau descrip- 
tion of the strong-coupling pairing, and the summary of the 
paper. 
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2. Model and Formulation 

We consider the Holstein-Hubbard model, 



H 



i/ep + 



= WE {b; + i j + g ^ [bi + bi) {nia - ^ j 



(1) 



ka 



where n,„ = c^I^c.^ is the electron density operator for the 
spin a, and the Einstein phonons b. and /?! with the frequency 
We couple with the electron density at the site /. As we men- 
tioned in the introduction, all the parameters are regarded as 
the renormalized ones. 

Using the path-integral form and integrating out the Ein- 
stein phonons, we obtain the partition function with the ef- 
fective action as 

a k 

S int = ^ X X ~ k')cjk)c/}(q - k)cp{q - k')Ca(k'), 

QjS kk' q 

(2) 

where k - {k, ia)„), 2* - Yjk,n- Here the effective interaction 
among electrons is given by 

-.2 



v(^) = tz-y,)- 



2g 

Vo = — , 
We 



(3) 



where - (q, e,„) with the bosonic Matsubara frequency e,„ - 
2nTm. 

Let us examine the general filling of the electron density n, 
and no particular structure of the Fermi surface is concerned 
for the superconductivity. In this case, the ordinary Cooper 
instability is expected to occur in the isotropic s-wave pair- 
ing. Following the recent theoretical development for the OF 
superconductivity, the gap equations for the i-wave EF 
singlet pairing t^apik) - A+(i(jj„){icr^')crfi and the i-wave OF 
unitary triplet pairing^"*' t^apik) - iA-(i(jj„)6ap are given by 

CO 

A0(/w„) = -r ^ V^{(x)„, u)n')F^ii(x)„'), ((p = +), (4) 

n'=0 

where we have used the symmetry A^(icj„) = (f)A^{-iojn) of 
the gap functions for the EF (((> = +) and the OF {(p - -) pair- 
ings {(A),, - nT{2n + 1) is the fermionic Matsubara frequency). 
Here we have introduced the interaction matrices as 

y±(w„, a)„>) - v{ia)„ - iun') ± v(/w„ + ito,,'). (5) 

The local anomalous green function is given by 

2pfA^(;w„) 



FJioj,,) = 



■ tan 



-I- |A^(/w„)p^ 



(6) 

Here we have used the constant density of states, p(^) - 
Pf S((^c - 1^1) for the conduction electrons and the cut-ofi" ujc 
is the order of the renormalized bandwidth. In order to deter- 
mine the transition temperature Tc, the linearized gap equa- 




Fig. 1. (Color online) The interaction matrices for (a) the EF channel, 
\V+(ai„,aj„,)\/Vo with U/Vq = 0.25, (b) the OF channel, \V-(a)„,cj„,)\/Vo 
in the upper panel. The lower panels show the diagonal parts, V+(ai„, oj„). 



tion is often used. 



AAJicon) = -2pFr V tan '1 — 1 AJia>„>), 

(7) 

where the eigenvalue A{T) has been introduced to check the 
development of the Cooper instability signaled by A{Tc) - 1. 

With the help of the solution of the gap equation (4), we can 
express the free energy^^* (measured from that of the normal 
state) as 

~ ( ^ojI + |A^(/w„)p - w„)^ 
T - -InppT y 



\A^{i(i)„) 



n=o ^ja)l + |A^(/w„)p ( + |A0(;w„)P + w,,) 

(8) 

The latter expression is convenient for the numerical compu- 
tation due to the fast convergence of the Matsubara summa- 
tion. The entropy and the specific heat can be computed by 
the numerical differentiation as 



dT ^ dS 
S^-^, C^T — . 
dT dT 

The superfluid density ^^-^^^ is given by 



ris - 2nT > , 



|A0(/w„)|2 



«=() + |A^(/w„)p) 



3/2' 



(9) 



(10) 



Hereafter, we use the unit of energy or temperature as 1 /2pF - 

Uc - I- 

3. Results 

Let us first clarify the difference of the interaction matrices 
between the EF and the OF channels. The upper panels of 
Fig. 1 show the intensity map of\V±{(jD„, (jj„')\IVq with U/Vq - 
0.25 for (a) the EF channel (+) and (b) the OF channel (-). 
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Fig. 2. (Color online) The T dependence of the eigenvalue of the linearized 
gap equation. For (/ = 0, (the open circle) for the EF pairing is slightly 
higher than that for the OF pairing. For = 0.08, no transitions occur in 
the OF channel. 



Fig. 3. (Color online) The ui,, dependence of the OF gap function. The inset 
shows the T dependence of A_(rai„=o). The dotted lines are a guide for 
eyes. 



For the EF channel, the dominant contribution comes from 
the square region of (jl>„, <y„' < and hence the summation 
in the gap equation (4) is usually cut off at oj„t ~ oje, and the 
diagonal part of lij„ - a)„' > oje is neglected. On the other 
hand, the diagonal part plays an essential role in the case of 
the OF pairing. Moreover, in the adiabatic limit we — > as 
consequences of strong fluctuations of ionic vibrations and/or 
a strong anharmonicity of the ionic potential, the diagonal part 
of the interaction matrices dominate the attractive interaction 
both for the EF and the OF pairings. 

The diagonal part of the interaction matrices given by 



y+(w„,w„) = 2f/-yo 



1 + 



1 



y_(w„,w„) 



1 



(2w„/we)2 + 1 



1 



(11) 



(2w„/we)2 + 1 

are shown in the lower panel of Fig. 1. Since the OF pairing 
can avoid the on-site Coulomb repulsion U, it is more sta- 
bilized against the EF pairing as a»E approaches to 0. At the 
same time, the strong-coupling condition igp = ppVo > 1 is 
necessary for an appearance of the OF pairing because the 
attraction is restricted to the diagonal part of the scattering 
(a>i,,—a>„) — > (a)„',—io„') with lj„' — co,, and we < w„ < 
|A_(/w„)| ~ TVq, yielding the net amount of the attractive 
contributions is rather weak. 

As T decreases, the dominant part of the attraction (oje < 
a>„ < \A-(i(L)„)\) for the OF pairing is vanishing. Hence, the 
OF pairing tends to exhibit a reentrant behavior^- '^^ On the 
contrary, the dominant part of the attraction for the EF pairing 
becomes stronger as T decreases, and it is possible to over- 
come the OF pairing in lower temperatures, provided the cou- 
pling constant is strong enough. 

Next, we discuss the T-oje phase diagram. For this purpose, 
we solve the linearized gap equation (7). Figure 2 shows the 
T dependence of the eigenvalues for g = 0.3. For loe - 0.05 
(/lep - 1.8) and U - 0, Tc for the EF pairing is slightly 
higher than that for the OF pairing. When U is switched on as 
U = 0.5, only the EF pairing is aff'ected to decrease T^, and 
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g=0.3, U=0.5, (i)£=0.065 
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Fig. 4. (Color onhne) The T dependence of the free energy for the OF (the 
solid line) and the EF (the dotted line) pairing. The veilical dotted line 
indicates the Ist-order transition temperature between the OF and the EF 
phases. 



the OF pairing emerges. As mentioned in the above, A(T) de- 
creases below T ~ a>E due to the vanishing attraction for the 
OF pairing. The reentrant lower Tc to the normal phase ap- 
pears in this case. For oje - 0.08 (/lep - 1.1), the eigenvalue 
does not reach unity and only the EF pairing occurs. 

The solution of the gap function for the OF pairing is shown 
in Fig. 3. Owing to the strongly retarded attraction, the gap 
function decreases very fast as w„ increases. The T depen- 
dence of the most dominant component, A^iinT), is shown in 
the inset of Fig. 3. It exhibits a maximum at T > we- Note that 
the fast decrease of the gap function in a)„ is also obtained in 
the case of the EF pairing, and it is the characteristic feature 
of the strongly retarded attraction. 

In order to determine a complete phase diagram, we should 
compare the free energy of the OF and the EF phases. The T 
dependence of the free energy (8) is shown in Fig. 4. Due to 
the reentrant behavior of the OF phase, there exists the Ist- 
order transition to the EF phase at T ^ 0.026 around which 
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0.05 0.06 0.07 0.08 0,09 



Fig. 5. (Color online) The T-lo^ phase diagram. The "1st" in the OF triplet 
phase indicates the 1 st-order transition to the EF singlet phase determined 
by the comparison of the free energies of both phases. The solid lines repre- 
sent the actual transition lines, while the dotted hnes indicate the fictitious 
transitions without the other phases. 
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Fig. 6. (Color onhne) The T dependence of the specific heat C{J)lyN 
where yfi is the T-linear coefficient in the normal phase. The disconti- 
nuity at Tc is much smaller than that of the BCS theory. The vertical line 
represents the Ist-order transition to the EF pairing. 



we expect a coexistent phase. 

By a similar analysis, we draw the uj-e-T phase diagram as 
shown in Fig. 5. The OF phase is more stabilized when the 
pairing interaction becomes strongly retarded with decrease 
of (We that simultaneously makes the coupling constant X^^ to 
increase. The upper Tc is mostly higher than we which reflects 
the strong-coupling condition necessary for the OF pairing. 
The reentrant behavior appears in a certain part of the phase 
diagram, although the EF pairing emerges as the low-T phase 
in most cases. Note that the EF phase should be suppressed 
as JJ increases. The OF phase exists only below iwe ~ 0.077 
(/lep S 1.17) indicating again that the strong coupling is nec- 
essary for the appearance of the OF pairing. 

The r dependence of the specific heat is shown in Fig. 6. 
The discontinuity ^C{Tc)IjnTc - 0.34 at Tc - 0.148 is much 
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Fig. 7. (Color online) The T dependence of the superfluid density, n^iT). 
Both for the EF and OF pairings, the slope of at Tc is rather small due 
to the characteristics of the attractive interaction. 

smaller than that of the BCS theory, 1 .43. The discontinuity 
for the EF pairing 0.24 at Tc - 0.042) is also rather small. 
Thus, the smallness of the discontinuity is paradoxically the 
consequence of the strong-coupling and the strongly retarded 
attraction. 

Figure 7 shows the T dependence of the superfluid density. 
The T dependence of n, near ns(T) ^ 1.18(rc -T)/Tc for 
the OF pairing and ^ 0.97(7^ - T)/Tc for the EF pairing are 
also smaller than that of the BCS theory, ns(T) ^ 2(Tc-T)/Tc. 
This is also due to the strongly retarded attraction. The rele- 
vance between the strongly retarded attraction and the small- 
ness of the specific-heat jump and the slope of the superfluid 
density will be discussed below. 

4. Discussions and Conclusions 

In the previous section, we have demonstrated the neces- 
sary condition for the appearance and the characteristic fea- 
tures of the OF pairing. The favorable conditions for the OF 
pairing are the strong-coupling and the strongly retarded at- 
tractive interaction. These conditions result in a Cooper pair 
whose extent is quite small in space and large in time. In this 
sense, it can be said that the OF superconductivity is the op- 
posite limit of the one in the original BCS theory. 

As was shown in the previous section, the w„=() component 
of the gap function plays a dominant role. Namely, for we ^ 
Tc, Ao (jJc and T < Tc, we adopt the Ginzburg-Landau 
description with the approximate gap function, 

A_(/w„) ^ /A()(5„,() - 6n-i), (12) 

(Ao is real and A_(t) - 2Aosin(7rrT) in the imaginary-time 
domain). In this approximation, the T dependence of Aq near 
Tc from the gap equation is given by 

A,iT)^nTc^^^C'l\ ^.lZ£.2(l-lj, (13) 

where - -ppV-{nTc,nTc) ~ /iep, and Aq > I. The quantity 
^ is regarded as the measure of the strong-coupling and the 
strongly retarded interaction. The free energy is then obtained 
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as 



2r2 



(14) 



which gives the discontinuity of the specific heat at T = as 



AC 



-r 



JnT, 

The superfluid density also reads 



JN 



2n pp 



= 2- 



(15) 



(16) 



These quantities exhibit the same exponents as those of the 

BCS theory, but the coeflicient is different by factors of cer- 
tain power of ^. It is interesting to note that the values of the 
BCS theory are almost reproduced by putting ( ~ I {Aq ~ 2). 
Thus, the coefficients become rather smaller when ^ < 1, 
while they become larger for ^ > 1 as in the case of the 
ordinary strong-coupling superconductivity. We should em- 
phasize however that these behaviors are the consequence of 
the strong retardation, and hence there exists little difference 
between the OF and the EF pairings near Tc, if the interac- 
tion is strongly retarded (Tc for the EF pairing is obtained by 
replacing V- in Aq with V+). A detailed analysis for the strong- 
coupling (adiabatic) limit of the Eliashberg theory (for the EF 
pairing) is available in the literature. ^''^ 

In this paper, we have neglected both the self-energy and 
the vertex corrections. Naively speaking, the self-energy cor- 
rection tends to decrease the effective coupling constant A^p 
by a factor of the mass renormalization, and the OF phase 
would be suppressed. As was shown in the previous section, 
the OF pairing is formed at rather high frequency (jj„=() = 
nT » a»E where the mass renormalization factor could be 
small as compared with the Iow-oj limit, 1 + A^^. Moreover, 
the strong retardation and the strong coupling are the neces- 
sary conditions for the OF pairing to appear In this situation, 
it has been argued the breakdown of the Migdal theorem^^^^' 
and the vertex correction should be treated properly. The in- 
clusion of the self-energy and the vertex corrections are left 
for the future investigation. 

In summary, we have discussed the possible OF super- 
conductivity in the strong-coupling electron-phonon systems 
based on the Holstein-Hubbard model. The presence of the 
particularly soft Einstein phonon mode provides the opposite 
situation of the original BCS theory, i.e., the strong-coupling 
and the strongly retarded attractive interaction. They are the 
favorable conditions for the OF superconductivity, where the 
extent of the Copper pair becomes small in space and large 
in time. The Ginzburg-Landau analysis for this situation con- 
cludes that the specific-heat jump and the slope of the super- 
fluid density at Tc are factorized by the ratio f = 4Tc/cOc, and 
they can be much smaller than those of the BCS theory for 
f <1. 
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